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Abstract 

Recently a renormalizable model of gravity has been proposed, which might be a UV comple- 
tion of General Relativity (GR) or its infra-red modification, probably with a strongly coupled 
scalar mode. Although the generic vacuum of the theory is anti-de Sitter one, particular limits 
of the theory allow for the Minkowski vacuum. In this limit (though without consideration of 
the strongly coupled scalar field) post-Newtonian coefficients of spherically symmetric solutions 
coincide with those of the General Relativity. Thus the deviations from the convenient GR should 
be tested beyond the post-Newtonian corrections, that is for a system with strong gravity at as- 
trophysical scales. In this letter we consider potentially observable properties of black holes in the 
deformed Horava-Lifshitz gravity with Minkowski vacuum: the gravitational lensing and quasi- 
normal modes. We have showed that the bending angle is seemingly smaller in the considered 
Horava-Lifshitz gravity than in GR. The quasinormal modes of black holes are longer lived and 
have larger real oscillation frequency in the Horava-Lifshitz gravity than in GR. These corrections 
should be observable in the near future experiments on lensing and by gravitational antennas, 
helping to constrain parameters of the Horava-Lifshitz gravity or to discard it. 

1 Introduction 

Recently P. Horava suggested a renormalizable four-dimensional theory of gravity, which admits the 
Lifshitz scale-invariance in time and space [I], [2]. Thus, except of a particular limit, the theory is with 
a broken Lorentz symmetry at short distances, while at large distances, higher derivative terms do 
not contribute and the theory runs to a kind of classical gravity theory, if a coupling, which controls 
the contribution of trace of the extrinsic curvature has a particular value (A = 1). An important 
observation that has been made recently is that in this infra-red limit of the theory, gravity needs to 
be strongly coupled to a scalar field [I], [3]. This happens because in the particular infra-red limit 
of the theory, the detailed balance is explicitly broken, instead of the full diffeomorphism the theory 
has a reduced one, and thereby should have more propagating degrees of freedom [7J than GR has. 
Various properties and consequences of the Horava gravities have been considered recently in a lot of 
papers [5]. 

In this letter we shall investigate a particular, non-minimal limit of the deformed Horava theory, 
which implies the Minkowski vacuum and a kind of modified gravity at infra-red, without consideration 
of the strongly coupled scalar mode. As in this specific limit, post-Newtonian corrections of the 
Horava Gravity coincide with those of the standard GR [6], the effects when the gravitational field 
is weak are well described by the post-Newtonian approximations. To test the theory in experiment, 
one needs therefore the astrophysically strong gravitational field effects. A system at hand is a 
supermassive black holes, whose parameters might be observed through its characteristic frequencies, 
called quasinormal modes [9]. The latter are proper oscillations of a black hole as a response to 
external perturbations. Another opportunity is to observe the gravitational lensing effects, such as 
bending angle and time delay, for a super-massive Galactic black hole at small impact parameters b. 



(At large b the lensing is in the regime of the weak field and thus well described by post-Newtonian 
approximations.) Here we shall consider both effects, quasinormal modes and gravitational lensing 
for a asymptotically flat black holes in the particular, non-minimal Horava gravity. We shall show 
that both effects give considerable deviation from the standard GR. 

The paper is organized as follows. Sec II gives the brief description of the black hole, under 
consideration. Sec III analyzes the geodesic motion for massless particles. Sec IV is devoted to the 
bending angle effect and Sec V gives quasinormal modes for the test scalar field in the background 
of the black hole. In the Conclusions we outline the number of open questions and summarize the 
obtained results. 



2 The black hole metric 

In the particular limit Aw — ► 0. the Horava-Lifshitz theory can be described by the following action 
S = J dtfzrfN {^{K^Sy - XK 2 ) - ^ Cij C« + ^e^VjR^ 

_^f R f) R m + 8( ^3 A /~/V 3) ) 2 + /^ (3) ) • (i) 



Here 



K l , = ^(g lJ -y i N J -V,N l ) , (2) 



(3) 



are the the second fundamental form and the Cotton tensor respectively. In the particular type of 
the Horava theory, we have 

c 2 = fcV/2, G = fe 2 /32vrc, A = 1. (4) 
The black hole solution for the A = 1 theory has the form (see formula (30) in [8] and [6]) 

ds 2 = -f{r)dt 2 + -— + r 2 (d9 2 + sin 2 ## 2 ), (5) 

where 

uj = 8/i 2 (3A - 1)/k 2 = uj = 16h 2 /k 2 (6) 



f{r) = 1 + ujr 2 - ^r(uj 2 r 3 + 4uM), (7) 



and M is an integration constant. 
There are two event horizons at 



r± = M [1 ± v 7 ! - l/2o;M 2 J . (8) 
Avoiding naked singularity at the origin implies 

uM 2 >\. (9) 

In the regime of the conventional General Relativity toM 2 ^> 1 the outer horizon approaches the 
Schwarzschild horizon r + ~ 2M and the inner horizon approaches the central singularity. 
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3 Geodesic motion 



In the Horava-Lifshitz gravity, there is no a full diffeo-invariance of the Hamiltonian formalism. In 
the theory with matter lagrangian an alternative covariant approach should be used instead [7], where 
the matter couples to gravity in a consistent way. Then corrected Hamiltonian equations should be 
used. Here, as a first step, we shall neglect these corrections and shall imply that at large distances 
and small momenta the ordinary Hamiltonian formalism can be used as some approximation. In this 
way, we include into consideration the influence of the astrophysical background (a black hole) on a 
light propagation, but do not include the deviation into the form of the propagator. The full analysis 
certainly should be based on the covariant approach to the Horava gravity [TJ, as well as on inclusion 
of strongly coupled scalar mode [3], [1]. This also concerns the quasinormal mode analysis related 
further. 

Our main aim will be analysis of the null geodesic motion with the help of the ordinary Hamiltonian 
formalism. However, we shall consider also some general features of both null and tim-like geodesies, 
implying therefore a massive term in the geodesic equation. The Hamiltonian has the form 

\ .... 1 o dx^ , . 

H = -gr Vp,Vu = -gA* > = 9v»-fo' ( 10 ) 

Then, the Hamilton-Jacoby equations for the time-like and null geodesies are 

1 dS_ dS_ = dS = _l , 

2 y dxudx" ds 2^ " { ' 

Here s is an invariant affine parameter. For the metric (5) this equation takes the form 

f ^fdS\ 2 f fdS\ 2 1 fdS\ 2 1 fdS\ 2 _ 2 

/ \dt) fydr) r 2 \de) rWflVW " M [ ' 

The action is: 

S = -i/i 2 s -Et + Lcj) + S r (r) + S e (6), 

and E and L are the particle's energy and angular momentum respectively, po = —E, and P3 = L, 
pi = S r and p2 = 5*6) are functions of r and respectively. Due-to spherical symmetry of the 
background any central plane of motion can be considered as an equatorial, and one can consider 
9 = tt/2 without loss of generality. 

Implying the normalization r = [is, where r is the proper time, the first integrals of motion are 

"^'l = E - "s = ± y< iJr ^' (13) 

The qualitative analysis of the motion can be made by considering the effective potential of the 
motion, 

U !ff = ^/(r) (l + ^) • (14) 

The effective potentials for various values of u are plotted on Fig. ([1]), where one can see that as u 
grows, the effective potential approaches the Schwarzschild form. 

It is well known that for the Schwarzschild black hole of mass M, there are two circular orbits of 
time-like geodesies, which radii are given by the equation dU e ff/dr = 0. The inner orbit corresponds 
to maximum of the effective potential and is therefore unstable, while the outer orbit corresponds to 
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a local minim and is stable. The minimal value of the momenta L, necessary for the existence of the 
outer stable orbit is 



or equivalently 



L > 2V3M/I. 



For null geodesies, there is no (stable) circular orbits in the Schwarzschild space-time, except for the 
one unstable orbit located in the maximum of the effective potential. 

The presence of oj slightly changes the situation. The radius of closed (unstable) circular orbit for 
null geodesies equals r c j r = 2.37228M for to = 1/2 and monotonically increases with u, approaching 
the pure Schwarzschild value r c j r = 3M for u>M 2 3> L (see Figj2]). Thus for null geodesies, as in the 
Schwarzschild case, there is only one unstable circular closed orbit at the maximum of the potential, 
though of smaller radius in the presence of u. 

In a similar fashion, the time-like geodesies around the Horava black hole has two circular orbits: 
an inner one (unstable) and an outer (stable). The condition for existence of the outer stable orbit 
reads 



L > 2 V / 3M^ 



where the function e(lo) is shown on Figj3j 



u. 

. 98 




Figure 1: The effective potential of the geodesic motion for time-like particles for L = 2y / 3M / u, 
M = fJ, = l,LJ = l/2 (red), 2 (green), and for the Schwarzschild black hole (blue). 



4 Bending angle 

The propagation equation for the null geodesies is 



.500(1 + ^2, h = L/E , 
911 \ 533 ' 



The geodesic trajectory equation for the null geodesies has the form 

V#/ 5ii V foo 



(15) 



(16) 



Thus the propagation and trajectory equations contain the ratio b, which is the impact parameter. 
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Figure 2: The radius of massless particle circular orbit r c j r as a function of co, M = 1. 
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Figure 3: The deviation from the Schwarzschild stability region of time-like geodesies is described by 
the function e(w)//iM. 
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When passing near the black hole, a ray of light approaches it at some minimal distance, called 
distance of closest approach r m i n . The latter can be found as the largest real root of the equation 
alr/dt = 0. For the asymptotically flat black holes (4-6) in the considered Horava-Lifshitz gravity, this 
equation reads 

r 2 + 6 2 ( v /rw(4M + wr 3 ) - r 2 to - 1) = 0. (17) 




Figure 4: The value of a + it as a function of u, M = 1, b = 6. 
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Figure 5: The value of a + ir as a function of u, M = 1, b = 10. 



Upon calculation of r m i n with good accuracy, one can perform integrations for finding bending 
angle: 

(18) 



a 



dr J r dr 



7T. 



Here r Q is radial coordinate of an observer and r s is radial coordinate of the source. One can see 
the bending angles as a function of uo on Fig. 3 and 4. for various values of the impact parameter. 
Starting from a conventional General Relativity, s super-massive black hole in the center of our galaxy 
has the mass M = (3.6 ± 0.2)10 6 Mq (what corresponds toMa 3.4 • 1CT 7 pc, in units c=G=l), and 
is situated at a distance r Q = (7.9 ± 0.4) 10 3 pc from the Earth. Typical distance from a source to 
a black hole is of order r s = 1 pc. Now let us take M = 1, then we can re-scale the corresponding 
values for r Q and r s 
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r Q = 2.3235 • 10 iU , r s = 2.9412 • 10 b , M = 1. 



(19) 



In the used here units, the bending angle is smaller for the Horava-Lifshitz gravity than for the 
convenient General Relativity. (However, in order to give some numerical estimations for experiments, 
one should include also the re-definitions of the fundamental constants (1)) In the regime of large 
ujM 2 the bending angle approach its Schwarzschild values. For relatively small u the deviation from 
Schwarzschild values is significant for small impact parameter b (see Fig. 3), while for large b, the 
deviations from Schwarzschild values are almost negligible. As large b regime is well described by 
the post-Newtonian coefficients of the metric, it confirms the fact that the metric (4-6) has the same 
Eddington-Robertson-Schiff Post-Newtonian parameters /3 = 7 = 1 as pure Schwarzschild metric 
does [10]. 

The observed here corrections would be measurable today if a pulsar were found to be lensed by 
the Galactic super-massive black hole. Also it might be detected with various planned experimental 
missions for observations of gravitational lensing. 



5 Quasinormal modes 

As in the case of lensing, our approach here suffers from a number of approximations, which we shall 
enumerate here: 1) We neglect the possible strong coupling of the scalar mode in the Horava gravity 
[1], [3]. We can do this as some approximation to the realistic picture because of the following reasons. 
Quasinormal spectrum of black holes consists of a set of modes with various multipole numbers t. 
In the regime of large I the quasinormal spectrum does not depend on the spin of the field under 
consideration and stipulated by the centrifugal term ~ f(r)£(£ + l)/r 2 in the effective potential. 
At the same time, for axial type of gravitational perturbations the background scalar mode is not 
excited. In other words, one can always select some modes of gravitation spectrum of a black hole, 
which will be well described by a test scalar field in the background of the black hole. This certainly 
implies that we neglected the influence of the scalar mode on the black hole metric itself. 2) We 
implied the existence of the approximate Hamiltonian formalism and corresponding consequences, 
such as existence of the traditional Green functions etc., because we expect that at large distances 
and relatively small momenta these corrections to the ordinary Lorentz-invariant physics should be 
relatively small. 

Thus the wave equations for the test scalar <3? is 

WV^g*^ = 0, (20) 

After separation of angular variable and introducing a new wave function \& s , the perturbation 
equation can be reduced to the wave-like form 

ill + (£2 _ Fs ( r ))^ = 0, dn = (A(r)B(r))-^ 2 dr, (21) 
where the scalar field effective potentials has the form: 



V 8 = A(r)^±^- + ±(*r)B-ir) + MrWir)). (22i 



Here A(r) = B{r) = f(r). The effective potential has the form of the positive definite potential 
barriers. The quasinormal modes are solutions of the wave equation (|2ip . with the specific boundary 
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Table 1: Fundamental (n = 0) WKB quasinormal frequencies for various values of to, compared to 



the Schwarzschild QNMs. 



u 


£ = 


£ = 1 


£ = 2 


0.5 


0.1213 - 0.0732 i 


0.3227 - 0.0723 i 


0.5343 - 0.0711 i 


1 


0.1208 - 0.0893 i 


0.3057 - 0.0891 i 


0.5048 - 0.0877 i 


1.5 


0.1171 - 0.0969 i 


0.3011 - 0.0925 i 


0.4971 - 0.0912 i 


2 


0.1151 - 0.0993 i 


0.2989 - 0.0940 i 


0.4935 - 0.0927 i 


3 


0.1134 - 0.1007 i 


0.2969 - 0.0953 i 


0.4901 - 0.0941 i 


5 


0.1129 - 0.1009 i 


0.2953 - 0.09632 i 


0.4874 - 0.0952 i 


10 


0.1113 - 0.1010 i 


0.2941 - 0.0970 i 


0.4855 - 0.0960 i 


Schwarzschild 


0.1105 - 0.1049 i 


0.2929 - 0.0977 i 


0.4836 - 0.0968 i 



conditions, requiring pure out-going waves at spatial infinity and pure in-coming waves on the event 
horizon. Thus no waves come from infinity or the event horizon, 

~ e ±ir *", r* -> ±oo. (23) 

Here we imply that u = Re(£u) + Hm(oj) , and ~ e lujt . Thus Im(ui) determines the damping rate 
of a quasinormal mode and Re(Co) determines its real oscillation frequency. There is no unstable 
(growing) modes for the test scalar field, as the effective potential is positive definite. Thus there is 
no modes with Im(uj) > 

It is well known that in many cases WKB method can give accurate values of the lowest (that is 
longer lived) quasinormal frequencies. The 6-th order WKB formula reads 

# = 0,1,2..., (24) 




where Q = cD 2 — V and the correction terms A, were obtained in [12], [13]. Here Q l means the i-th 
derivative of Q at its maximum with respect to the tortoise coordinate r*. This formula was effectively 
used in a lot of papers (see for instance [14] and references therein). Fortunately, for the considered 
here case, WKB series shows convergence in all sixth orders. In general, this is not guaranteed because 
the WKB series converges only asymptotically. 

On table I and II one can see that for relatively small u the Re(ui) is larger than in the Schwarzschild 
limit, while the damping rate —Im{Cj) is smaller than for pure Schwarzschild. Thus in the considered 
Horava theory, quasinormal modes are longer lived and have larger oscillation frequency. In other 
words, a black hole in the particular Horava theory is a better oscillator than a Schwarzschild black 
hole. In the limit ujM 2 3> 1, quasinormal modes approach their Schwarzschild values (Table I, II). 
The approaching of the Schwarzschild limit can be better seen for the fundamental overtone, because 
the WKB approximation works better when £ > n, though, on Table II one can see approaching the 
Schwarzschild limit for the relatively small overtone numbers n. The pure Schwarzschild QNMs given 
in the tables I and II are accurate, because they were computed by the convergent Frobenius method. 



6 Conclusions 

In the present paper we have considered the quasinormal modes and the bending angle around non- 
minimal asymptotically flat black holes in the Horava gravity. We have shown that 
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Table 2: Higher overtones, calculated by WKB formula, for various values of uj, compared to the 
Schwarzschild QNMs, 1 = 2. 



n 


uj = 1/2 


w = 1 


uj = 3 


Schwarzschild 





0.5343 - 0.0711 i 


0.5048 - 0.0877 i 


0.4901 - 0.0941 i 


0.4836 - 0.0968 i 


1 


0.4748 - 0.3625 i 


0.4892 - 0.2665 i 


0.4720 - 0.2871 i 


0.4639 - 0.2956 i 


2 


0.5146 - 0.2147 i 


0.4617 - 0.4540 i 


0.4411 - 0.4925 i 


0.4305 - 0.5086 i 


3 


0.4146 - 0.5184 i 


0.4281 - 0.6518 i 


0.4061 - 0.7136 i 


0.3939 - 0.7381 i 


4 


0.3346 - 0.6872 i 


0.3943 - 0.8596 i 


0.3762 - 0.9493 i 


0.3613 - 0.9800 i 



1) The quasinormal modes of the non-minimal asymptotically flat black holes in the Horava gravity 
are longer lived and have larger real oscillation frequencies than their Schwarzschild analogs. 

2) The bending angle is seemingly smaller in the considered Horava gravity than in the standard 
GR and the deviation might be detectable if a pulsar were found to be lensed by the Galactic super- 
massive black hole or with the planned missions. 

In the near future experiments, this may help in constraining of the coupling constants in the 
Horava gravity. 

Let us enumerate here a number of ways, in which our work may be improved. 

First of all, in a similar fashion one can find the time delay, which is the difference between the 
light travel time for the actual ray, and the travel time for the ray the light would have taken in the 
Minkowskian space-time: 

rr min dt rr df rf 

t s -t = - Tdr+ dr-^-. (25) 

L ar J, r . dr cos is 

" ' 3 u I mm 

Here the term represents the propagation time for a ray of light, if the black hole is absent 
(see for instance |llj). However the direct integration may be unstable for the considered case, and 
a transformation of the integral may be necessary. One can also generalize the geodesic analysis to 
massive particles, as well as to particles with spin. 

The search for quasinormal modes can be continued to fields with non-zero spin, and, what is 
more important, to gravitational perturbations. The latter will test the stability of black holes in the 
Horava theory and is therefore an appealing problem. The scalar quasinormal modes, considered in 
this paper, should give qualitative description of the gravitational ones, at least in the regime of high 
real oscillation frequencies (eikonal regime), that as we already mentioned does not depend on a spin 
of a field. 

Note added. On the same day, when the first version of this manuscript was published in arXiv, a 
complementary paper of Songbai Chen and Jiliang Jing 15J on quasinormal modes of the deformed 
Horava-Lifshitz gravity appeared in arXiv. 
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